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Abstract 

In this fluid dynamics video prepared for the APS-DFD Gallery of 
Fluid Motion we study the Kelvin-Helmholtz instability for both New- 
tonian and viscoelastic jets. The nonlinear dynamics of the jet motion 
are slowed down by orders of magnitude using a synchronized strobe ef- 
fect coupled with precise timing control of perturbation frequencies. Our 
results show that at high wave-numbers the imposed perturbations ini- 
tially grow linearly with time and the jet axis remains straight while the 
Kelvin-Helmholtz wave amplitude grows and rolls up into bags that en- 
capsulate the central jet within themselves. At low wave-numbers (long 
wave-lengths) the jet axis buckles under the action of viscous stresses 
and a coupling between the Kelvin-Helmholtz instability and bending of 
the jet leads to new concertina or chevron modes which grow with time 
as they move downstream. Addition of viscoelasticity to the jet leads 
to the pronounced inhibition of the Kelvin-Helmholtz instability as the 
jet perturbation amplitude grows and large elongational stresses in the 
fluid become important. For long waves, the initially-relaxed viscoelastic 
jet first buckles in a manner similar to the Newtonian solvent but once 
again the viscoelastic effects suppress the instability growth as they are 
convected downstream. 

1 Introduction 

Kelvin-Helmholtz instability is a well-known phenomenon, which occurs in na- 
ture when two identical or dissimilar liquids with different velocities meet each 
other at an interface. Wave generation on water surfaces in a windy day, shear 
layer patterns in clouds, and even wave patterns observed in solar winds repre- 
sent common examples showing this instability in nature [1]. Kelvin-Helmholtz 
instability also occurs in the industrially-important process of air-assisted atom- 
ization (between two coaxial jets) as the primary instability generating waves on 



the liquid surface [3]. Helmholtz (1868) was the first scientist who noted this in- 
stability and remarked: every perfect geometrically sharp edge by which a fluid 
flows must tear it asunder and establish a surface of separation, however slowly 
the rest of the fluid may move [1]. Lord Kelvin (1871) later developed this idea 
furthermore and established a mathematical solution for the problem. He later 
used the results to understand wave generation by wind blowing across a flat 
ocean surface. The effects of viscoelasticity on the Kelvin-Helmholtz instability 
are still largely unknown, despite the fact that viscoelastic jets and interfaces are 
key components of many industrial applications such as air-assisted atomization. 
The current study investigates the effects of viscoelasticity on this instability, in 
an axisymmetric flow configuration, benefitting from a novel strobe-based high 
resolution experimental setup. 

2 Experimental Setup 

The setup used in these experiments (Figure 1) has a sub- millimeter nozzle 
(diameter: 150 /am ) from which fluid is discharged into a water tank. The tank 
is made of transparent flat walls, which provide optical access to the axial po- 
sition where instability starts. The nozzle is perturbed using an annular piezo- 
electric actuator squeezing the nozzle walls within a wide range of frequencies 
(1Hz < f < <10kHz ). The nozzle and the tank are both positioned on an two- 
axis positioning stage, mounted on a stepper motor for precise motion in the 
axial direction (Figure 1). Liquids are pumped using a syringe pump (4400 Ultra 
PHD from Harvard Apparatus) at a fixed rate with precise control on the pump- 
ing rate (errors in flow rate are less than 0.5 percent). Two signal generators are 
used for perturbing the jet; one carries the signal to the piezoelectric actuator at 
moderately small amplitudes (1-5 volts) and the other generator sends a similar 
signal with slightly different frequency / — Sf to the strobe light. This slight 
difference in frequencies between the perturbation and the pulsed lighting leads 
to an apparent slowing down of the motion by a factor of Sf/f . 

3 Test Fluids 

The background liquid inside the tank is DI water, which is kept at room tem- 
perature. Two different liquids are discharged in the water tank; a Newtonian 
solvent (water and glycerol mixture) and a dilute viscoelastic solution made by 
dissolving a small amount (c = O.lpercent) of Poly-Ethylene-Oxide (PEO-300K) 
into the same Newtonian solvent. The shear viscosities of these two liquids are 
very close to each other at high shear rates (a common property of dilute PEO 
solutions) also their liquid/air surface tensions are almost identical since PEO is 
not active at the surface. The only difference between these two is the addition 
of extensional viscoelasticity, which results from dissolving PEO polymer chains 
in the solvent [4]. The material parameter that quantifies this difference is the 
extensional relaxation time (A# ). For the Newtonian solvent tested (like any 
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Figure 1: The experimental set-up used in this study. 



other Newtonian liquid) the relaxation time is zero since there is no intermolec- 
ular mechanism for elasticity or stress relaxation, however by addition of 0.1 
percent of PEO the relaxation time increases to A# = 0.36 ms (measurements 
were done using jet break-up extensional rheometry [5]). 

4 Contents of the Video 

As the Newtonian liquid leaves the nozzle and discharges into the water tank 
with a velocity around 3 m/s, shear instabilities start to appear due to the 
difference in velocities at the liquid-liquid interface. The amplitude of the wave 
crests grow with time as the jet moves downstream and become elongated in the 
upstream direction as they protrude into the stationary outer fluid; then due to 
the axisymmetry of the geometry, the wave crests form re-entrant bags that wrap 
around the jet axis entraining the outer fluid. The captured movies show that 
at high wave-numbers (or short wavelengths) the Kelvin-Helmholtz instability 
develops for the Newtonian solvent right after the nozzle exit and the instability 
evolves rapidly into a nonlinear state consisting of stacked chevron-like waves. 
At these high wave-numbers the jet remains axially straight throughout the 
evolution of the Kelvin-Helmholtz instability and there is no sign of bending or 
buckling in the jet axis (Figure 2). At low wave-numbers (long wavelengths) for 
the Newtonian liquid the shear instabilities start again close to the nozzle exit 



3 




Figure 2: Three snapshots of the Kelvin-Helmholtz instability (with 40 ps time dif- 
ference) for the Newtonian liquid and short perturbation wavelengths. 

kRo = (2tt//V) Ro = 0.8 , p /p w = 1.103 , p/p w = 3.2 , 
Re = (pVD/p) = 146 , Oh = p/Vp^D = 0.1 , Wi = X E V/D = 



but at the same time the jets axis starts to buckle due to the viscous compressive 
stresses acting on the fluid column [6]. Bending of the jet, coupled with the 
nonlinear evolution of the waves resulting from Kelvin-Helmholtz instability, 
leads to a complex sinuous wave pattern of growing amplitude (which has the 
same frequency as the imposed perturbation). The buckling and associated 
deceleration in the viscous jet concertinas each wave crest into the next one 
in a shape very similar to what one observes in tulip petals. These tulip like 
patterns are convected downstream and grow approximately linearly in radius 
with time while becoming progressively compressed and stacked together in the 
z-direction (Figure 3). 

For the viscoelastic jet, at high wave- numbers (short wavelengths), the Kelvin- 
Helmholtz instability again starts from the nozzle and in the linear region the 
instability behaves very similarly to the Newtonian solvent. However, as the 
instability grows to larger amplitudes the wave crests show enhanced resistance 
to the elongation due to the presence of viscoelasticity and this consequently 
inhibits the nonlinear bag formation and encapsulation process (Figure 4). At 
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Figure 3: Three snapshots of the Kelvin-Helmholtz instability (with 40 ps time dif- 
ference) for the Newtonian liquid and long perturbation wavelengths. 

kRo = (2tt// V) Ro = 0.3 , p/p w = 1.103 , p/p w = 3.2 , 
Re = (pVD/p) = 146 , Oh = p/Vp^D = 0.1 , Wi = X E V/D = 
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Figure 4: Three snapshots of the short-wavelengthKelvin-Helmholtz instability (with 
40/xs time difference) for the viscoealstic liquid. 

kRo = (2tt//V) R = 0.8 , p /p w = 1.103 , fi/fi w = 3.2 , 
Re = (pVD/p) = 146 , Oh = p/Vp^D = 0.1 , Wi = X E V/D = 



short wave- numbers (long wavelengths) the viscoelastic jet again undergoes a 
buckling instability coupled with Kelvin-Helmholtz perturbations very similar 
to the Newtonian jet (Figure 5) but the tulip patterns are not convected down- 
stream in a regular way and the width of the jet grows much more slowly than 
was observed in the Newtonian solvent due to the additional viscoelastic stresses. 
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Figure 5: Three snapshots of the Kelvin-Helmholtz instability (with 40 ps time dif- 
ference) for the viscoelastic liquid and long perturbation wavelengths. 

kRo = (2tt//V) Ro = 0.3 , p/p w = 1.103 , p/p w = 3.2 , 
Re = (pVD/p) = 146 , Oh = p/Vp^D = 0.1 , Wi = X E V/D = 
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